Optimal Prediction for Moment Models: Crescendo Diffusion and Reordered Equations

A direct numerical solution of the radiative transfer equation or any kinetic equation is typically expensive, since the radiative intensity depends on time, space and direction. An expansion in the direction variables yields an equivalent system of infinitely many moments. A fundamental problem is how to truncate the system. Various closures have been presented in the literature. We want to study moment closure generally within the framework of optimal prediction, a strategy to approximate the mean solution of a large system by a smaller system, for radiation moment systems. We apply this strategy to radiative transfer and show that several closures can be re-derived within this framework, e.g. PN , diffusion, and diffusion correction closures. In addition, the formalism gives rise to new parabolic systems, the reordered PN equations, that are similar to the simplified PN equations. Furthermore, we propose a modification to existing closures. Although simple and with no extra cost, this newly derived crescendo diffusion yields better approximations in numerical tests.